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Abstract 

We discuss the possibility of representing supersymmetry exactly in a 
lattice discretized system. In particular, we construct a perfect supersym- 
metric action for the Wess-Zumino model. 

1 Introduction 

Lattice simulations of supersymmetric systems usually apply formulations, which 
reveal the supersymmetry in the continuum limit but not in the lattice discretized 
version IJ. This note addresses the issue of a lattice action, which directly displays 
a continuous form of supersymmetry. Some related works are listed in Ref. 0. 

1.1 A simple supersymmetric model 

We first consider a simple 2d SUSY model |J given in the continuous Euclidean 
space by the Lagrangian 

C = ^j^ip + d^d^if , (1.1) 

where ip is a real "Majorana spinor" [] and <p is a real scalar field. Many qualitative 
features with respect to a lattice formulation are the same as in the Wess-Zumino 
model. The action S is invariant under the field transformations 

Sip = -y^p e , 5ip = eip, (1.2) 

1 Strictly speaking, there are no Majorana spinors in Euclidean space, we just refer to the Eu- 
clidean version of the corresponding formulae in Minkowski space. A definition is given for instance 
in Ref. j|, based on a Euclidean analog of charge conjugation. Note that 4> and ip are not 
independent. 



where the two components of the transformation parameter e are real Grassmann 
variables. As an important property, we note that the supersymmetric generator 
forms a closed algebra with the translation operator, 

[6 1 ,6 2 )<p = (ei 7 M e 2 - e 2 7^ ei) . (1.3) 
1.2 Ansatz for a lattice formulation 

Let us consider a rather general ansatz for a lattice discretization of this free system 
in momentum space, 
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s = Lwp {*(-p)[^(p) + Kp)Mp) + H-pMpMp)} , 

8*(p) = -[ 7M i2 M (p)+L(p)]$(p)e, 

= e[«(p)+7^(p)]* . (1.4) 

Here, ^ and $ are the massless lattice fermion resp. scalar field, and the new 
quantities, which we introduce as an ansatz for the inverse propagators and for 
the transformation terms A, tt, L, u, v^) are real in coordinate space. It is 
desirable that they are all local, i.e. analytic in momentum space. We require the 
low energy expansion of the action to reproduce the correct continuum limit, and 
the inverse propagators obey in coordinate space the normalization conditions 

5> M p„(z) = 1 , £A(x) = , 5> 2 ft(;r) = -4 . (1.5) 

X XX 

We do not require the lattice transformation terms to correspond exactly to the 
continuum transformations (|1.2|) . Hence this ansatz includes a possible "remnant 
supersymmetry" of the lattice action, similar to the Ginsparg- Wilson relation for 
the chiral symmetry |J. The general (remnant) supersymmetry requirement SS = 
amounts to 

- Pfl (-p)] + L(-p)[X(-p) - A(p)] + u(p)[Q(p) + Q(-p)] = , 
R^-pM-p) - A(p)] + L(-p)ip^-p) - Pli (p)\ + v^p)[Q(p) + Q(-p)] = . (1.6) 

It is sensible to assume the following symmetry properties in the action: p M is odd 
in the \i component and even in all other directions, while the Dirac scalars A and 
Q are entirely even. Then the conditions simplify to 

-R^-p) P ^p) + u(p)n(p) = 0, 

-L(-p)pr(p)+Vr(p)n(p) = 0. (1.7) 
Remarkably, A does not occur any more in these conditions. 



2 In the first expression, there is no sum over /i. 
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Finally, we assume for the transformation terms and L the same symmetries 
as for p M , A, respectively, which leads to 



- R»{p)Ph{p) = u{p)fl(p) , 

L(P)PM = v,(pMp) = 1,2). (1.8) 

The translation operator is identified from 

y 

Q(x-y) = Y,\. u ( x ~ z ) + lu"0u{x - z)\ [y^R^z - y) + L(z - y)\ , (1.9) 

z 

but this general form is not immediately instructive. 

Let us consider simple solutions for the case u = 1, = L = 0. f] The standard 
lattice action, p^{j>) = ip^ : = isinp M , Vt{p) = p 2 := ^^[2 sin(p Ai /2)] 2 , requires 

p 2 

Rpip) = -=o iPn > 

which is singular at p — (vr, 0), (0,7r) and (tt,tt), hence the transformation is non- 
local. An obvious concept to simplify - and to obtain the same dispersion relation 
for fermion and scalar - is the use of the same lattice differential operator for the 
scalar and the fermion part of the action. 

One way to do so is to set fl(p) = p 2 , R^ip) = P^p) = ip^, which is local but 
affected by doubling, both, for the fermion as well as the scalar. In the present 
model, unlike the case of the Wess-Zumino model ||, we cannot treat them by 
adding Wilson terms ((r/2) -p 2 ), because terms of this kind alter Q but not p M (the 
fermionic Wilson terms contributes to A). Hence R^ gets complicated and non-local 
again, R^{p) = ip^{l + (r/2)p 2 /p 2 ). 

If one is ready to accept non-locality, then it looks simpler to adjust the differential 
operators the other way round, p^{p) = R^ip) = ip^, as suggested in Ref. 0, and 
£l(p) = p 2 . Then the translation operator resulting from the lattice version of eq. 
( |1.3| ) corresponds to a half-lattice shift, whereas it is a full lattice shift for the option 
mentioned before. However, the fermionic inverse propagator performs a finite gap 
at the edge of the Brillouin zone, so we are dealing with a non-locality similar to 
the SLAC fermions. This suggests that also this approach fails to recover Lorentz 
invariance in the presence of a gauge interaction, as was pointed out for the SLAC 
fermions on the one loop level ||. 

One can construct a number of solutions of this type by hand. For instance, the 
standard action together with u(p) = Yl^ cosQo^/2) even provides locality, but such 
hand-made constructions look hardly satisfactory. Similar to the problem of chiral 
fermions on the lattice, they do not appear promising for the consistent incorporation 
of interactions. Hence we are going to follow a different, more systematic strategy. 

3 Note that the two transformation terms with an unusual Dirac structure, L and can only 
come into play simultaneously. 
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2 A perfect supersymmetric lattice action 



Since we are considering a free theory here, we can construct a perfect lattice action 
by "blocking from the continuum" , which corresponds to a block variable renormal- 
ization group transformation (RGT) with blocking factor infinity, 



e -3L*,»j = I Dl p D(p e -,wm-T\v,i>,*M , (2.1) 

where S is the perfect lattice action (i.e. an action without lattice artifacts), s the 
continuum action, and T the transformation term. We choose the latter such that 
the functional integral remains Gaussian, 

T = EI**" / ip(u)du} (af)^ [% - [ ^{u)du} 

x,y J C X JC y 

+ 5>*-/ <fi(u)du] (a% [$ y - [ <p(u)du}, (2.2) 



where C x is the unit square with center x, and cr , a s are arbitrary RGT parameters 
(a s has to be positive). The resulting perfect action reads 



\2 



n(» + 27r/) 2 , AS/ . ^Il(» + 27r/) z 

A ' (p) = S i^fk + Q/(p) • = g w + Q (p) • (2 - 3) 

where II (p) := Yl^p^/p^- Locality requires a* ^ 0, which naively breaks the chiral 
symmetry. However, the latter is still present in the observables ||, and a continuous 
remnant form of it even persists in the lattice action, if cr (p) is analytic ||. [] 

Now we consider the SUSY transformation. The variation of the continuum fields 
is given in eq. ( |1.2|) . If we transform simultaneously the lattice fields as 



5V X = -7 M / d^(u)du e , 8<$> x = e / i/;(u)du , (2.4) 

J (J x J (J x 

then all the square brackets in the expression for T (eq. ( |2.2j )) remain invariant - 
and so does the continuum action - hence SS = 0. 

Everything is consistent since we block the fields as well as their variations from 
the continuum. Note that this is not a solution along the lines of section L2, because 
the transformations of the lattice fields are not expressed directly in terms of lattice 
fields. (The special case of a 5 function RGT, a s , a? — > 0, is an exception, see 
below.) Hence the solution is somehow implicit. 

However, we can re-write the field variations in terms of lattice variables. First, 
we define a continuum current 

3ii = ■ (2-5) 

4 In this case, the full chiral symmetry of the fermion propagator is broken only locally. This is 
the property denoted as Ginsparg- Wilson relation. 
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We now block this current by integrating the flux through the face between two 
adjacent lattice cells, 



f 1/2 1 

Jn,x = / + + U ») du » \ v ^ A*» A = !)■ ( 2 - 6 ) 



This is a perfect lattice current 0. Here we assume it to be implemented so that 
eq. ( |2.6| ) holds exactly. As an interesting property, its lattice divergence is equal to 
the blocked continuum divergence, 

V^J^ := XX^,z ~ ^avs-a) = / d nU( u ) du (Gauss' law). (2.7) 

We are now prepared to write the variation of the lattice fermion field from eq. Q2.4Q 
in terms of lattice variables, 

5tf = -V M J M e , = -e V M J M . (2.8) 

In addition, we introduce the fermionic lattice field 



# x := / ij(u)du , (2.9) 

J C x 

which allows us to write also 5$ in terms of lattice quantities, 

5$ = e^. (2.10) 

For a 5 function RGT in the fermionic sector we have — ^f, but for finite this 
does not hold exactly. 

A generalization is possible for instance with respect to the blocking scheme. 
Instead of the block average scheme we have used so far, we can start from a more 
general ansatz 

T = - [li f (x- u)i){u)du] (a f )~^ [V v - [u f (x - u)^(u)du] 

x,y J J 

+ Et $ -- / n s (x-uMu)du] (a*)" 1 [%- f U s (x-uMu)du] , (2.11) 

where JU^(u)du = JU s (u)du = 1. Both convolution functions IF, II s are peaked 
around zero and decay fast enough to preserve locality on the lattice. 
Correspondingly, the variations of the lattice fields turn into 

5^ x = -7 M J TL f (x - u)d^(u)du e , 

5$ x = eju s (x-u)ip(u)du:=e^ x , (2.12) 

where we have adjusted the definition of \F If we want to achieve ^ — ^f, then we 
need - except for w — - also IF = IF. 
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This generalized scheme does not yield an obvious lattice current any more; the 
latter is a virtue of the block average scheme (characterized by IF (u), IP (u) = 1 if 
u 6 Co, and otherwise). In the general case, it is easier to consider the continuum 
divergence 

d(u) = dfj^u) = "fn 9 M u ) > ( 2 - 13 ) 
and build from it directly the lattice divergence 

D x = J IF" \x - u)d{u) du , 5^/ = -De. (2.14) 

Regarding the transformation algebra, we obtain 

[Si,S 2 ] $ 2 = (ei 7 M e 2 - e 2 7^ e x ) / U s (x - u)d^(u) du . (2.15) 

In particular, for the block average scheme this simplifies to 

[Si, 5 2 } $ = ei V /t J M e 2 - e 2 V M J M e a . (2.16) 

We see that the continuum translation operator is inherited by the perfect lattice 
formulation in a consistent way: eqs. ( |2.15| ), (|2.16|) show that the corresponding 
lattice translation operator is just the blocked continuum translation operator. The 
formula for [<5i,<5 2 ]\I/ is analogous. If we require the resulting translation operators 
to be identical, then we need IF = II s . Then the algebra of field variations and 
translation closes under the blocking integral. 

In any case, the fermionic and scalar spectrum are equal, because they both co- 
incide with the continuum spectrum. 

It is straightforward to apply this perfect lattice formulation to more compli- 
cated cases, see below. Interactions can be included perturbatively in the process of 
blocking from the continuum. For asymptotically free theories, at m = even the 
classically perfect action behaves perfectly [fLOR . Hence by means of an implicit (but 
not just symbolic) definition of the action - in terms of classical inverse blocking - 
we can also go beyond perturbation theory in the massless case. This is analogous 
to the fixed point formulation of a chiral gauge theory on the lattice P, |j~T|. 

3 Adding an auxiliary scalar field 

To proceed to the 2d Wess-Zumino model, we include an auxiliary scalar field.f] This 
equilibrates the number of fermionic and bosonic degrees of freedom. The continuum 
Lagrangian and the field variations read 

£ = ^c^V + d^ipdf.cp + f , 

SiP = -{j^d^ip + f) e , 5ip = e ip , 5f = e i^d^ip . (3.1) 

5 This is even necessary for the general validity of the resulting translation operator; otherwise 
it only holds on-shell. See for instance P. Freud, "Introduction to Supersymmetry" , Cambridge 
University Press, 1986. 
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Now the method presented in Ref. || is applicable in an extended form, if we 
use the following lattice discretization: 



i2 

S = Jb (2^ [*(-p)*7mP/i*(p) + ®(-p)P 2 ®(p) + F(-p)F(p) 

+ty(-p)W f (p)ty(p) + 2®(-p)W s (p)F(p) - $(-p)W a (p) 2 $(pj\ , 
6*(j>) = -{^P, + W s (p)Mp) + F(p)}e 
5$(jp) = e*(p) 

5F(p) = e[i ltl p^-W s {jp))^{jp) , (3.2) 

where W*(p), W s (p) are some sort of Wilson terms (zero at the origin, non-zero at 
the edges of the Brillouin zone, local and 2ir periodic, which implies that they are 
even). Hence they remove the degeneracy of the physical particles with their dou- 
blers. The standard form 1/2 p 2 is an example, but we can also insert more general 
scalar and fermionic Wilson terms [] and we always arrive at 5S = 0. If we also want 
the fermion and scalar spectrum to coincide, then we have to relate and W s . 
The procedure applied in Ref. || further restores a remnant chiral symmetry by 
means of the so-called overlap formalism, and this could also be done here. 

Instead we can apply the perfect action machine from section ^|, starting from 
the continuum system ( |3.1| ). This also solves the doubling problem and maintains a 
(remnant) chiral symmetry. The perfect propagator of the lattice field F reads 

A s (p) = £ IT% + 2ttZ) 2 + a\p) , (3.3) 
zez 2 

where a s is a RGT parameter analogous to a s . We now introduce two continuum 
currents, 7^y? and 7 M ^>, and we construct perfect currents from them. The explicit 
formulae are a straightforward extension of the formulae in section |[ We do not 
display them here, but we write them down for the further extension to the 4d 
Wess-Zumino model in the appendix. 



4 Conclusions 

We illuminated the problem of a direct construction of a supersymmetric lattice 
formulation. Then we have shown how a construction in terms of renormalization 
group transformations can be achieved. We preserve invariance under a continuous 
supersymmetric type of field transformations in a local perfect lattice action, which 
has also a remnant chiral symmetry. This applies to the 2d models discussed above, 
as well as to the 4d Wess-Zumino model, see appendix. We remark that the perfect 



6 The continuum limit of this action does in fact correspond to eq. (3.1), as we see if we substitute 
F{p) — > F(p) = F(p) + §(p)W s (p). Then the mixed term of F and <I> disappears, and we obtain 
another irrelevant term — $(— p)W s (p) 2 $(p). 
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formulation also cures the well-known problems related to the Leibniz rule ' - 
which breaks down for usual lattice discretizations - because here we keep track of the 
exact continuum differential operators. This is manifest in the translation operator, 
which results from a commutator of field variations: in the perfect lattice formulation, 
we obtain the consistently blocked continuum translation operator. Therefore the 
algebra with the field variations closes. 

Moreover, in the perfect lattice formulation, the fermionic and scalar dispersion 
relation coincide automatically. 

The next step is the inclusion of the gauge interaction; this work is in progress. A 
consistent blocking of the gauge fields from the continuum leads to a perfect action 
with all the continuum symmetry properties in the observables - and also in the 
action, if the transformation term respects these symmetries - but this construc- 
tion can only be performed perturbatively. For asymptotically free theories in the 
massless classically perfect action - constructed by simplified inverse blocking 

(based on minimization) - is sufficient for the same purpose and enables the step 
beyond perturbation theory. 

Acknowledgment I am indebted to M. Peardon for many helpful discussions. I 
also thank him and I. Montvay for reading the manuscript. 

A Application to the 4d Wess-Zumino model 

For completeness, we write down the corresponding formulae for a perfect lattice 
formulation of the 4d Wess-Zumino model: 



J Dil>D<pV>D<pWDfWDfW e 



,$(2)^(1)^(2)] _ / n J,n M (l) nm (2) nf (l) nf (2) - S [^(lV2),/(D,/(2>] 



X 



exp{-T[*, V>, $«, y,« $< 2 \ <p®, F {1 \ f^}} , (A.l) 

<#, ^\^\ / (2) ] = ~ frx [faM + (d^r + (^ (2) ) 2 + 2 + / (2) 2 ] , 

x,y C x C y 

2 

+ - / V {i) {u)du] (at)~i - / <p®(u)du] 

i=l x,y J C X JC y 

+[Ff - / f { Hu)du] (at)~; [FW - / /«(„)<*„] } , (A.2) 



' C x 

where a? and a\ have to be positive, and we are using the block average scheme. 



The perfect action S can be assembled by analogy from eqs. (|2.3| ) and (|3.3| ). The 
transformation of the continuum fields have the usual form, 

W = -[7m(^ (1) +75^ (2) ) + / (1) + 75 / (2) ]6, 
5(f^ = e if) , 8tp^ = e 75^ , 



7 See in particular the first paper by S. Nojiri. 



8 



5f 



(i) 



and the lattice field transformations amount to 
6V X = - 



6F® 



(A.3) 

{^^(u) + l5 d^ 2 \u)] + /«(«) + l5 f {2 \u)}du e 

V,JJ^ + JjJ + + 75^ (2) ] e, 
^(m)c?m = e , 5$i 2) = e 75 / ^(«)dw = e 7s^x , 

e 7^ / d,j,ip(u)du := e V^i^, 

e 757m / d^(u)du = e 7 5 V jU J (U , a; , (A.4) 



where J M , and are perfect lattice currents. 

Again the resulting translation operator is the blocked continuum translation 
operators, it forms a closed algebra with the field variations, the Leibniz rule is 
satisfied, and the dispersion relations of \I/ and $ coincide. 
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